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ABSTRACT 

We present new spectroscopic observations that are part of our continuing monitoring campaign of 88 quasars 
at z < 0.7 whose broad H/3 lines are offset from their systemic redshifts by a few thousand km s“*. These 
quasars have been considered candidates for hosting supermassive black hole binaries (SBHBs) by analogy 
with single-lined spectroscopic binary stars. We present the data and describe our improved analysis tech¬ 
niques, which include an extensive evaluation of uncertainties. We also present a variety of measurements 
from the spectra that are of general interest and will be useful in later stages of our analysis. Additionally, we 
take this opportunity to study the variability of the optical continuum and integrated flux of the broad H/3 line. 

We compare the variability properties of the SBHB candidates to those of a sample of typical quasars with 
similar redshifts and luminosities observed multiple times during the Sloan Digital Sky Survey. We find that 
the variability properties of the two samples are similar (variability amplitudes of 10-30% on time scales of ap¬ 
proximately 1-7 years) and that their structure functions can be described by a common model with parameters 
characteristic of typical quasars. These results suggest that the broad-line regions of SBHB candidates have a 
similar extent as those of typical quasars. We discuss the implications of this result for the SBHB scenario and 
ensuing constraints on the orbital parameters. 


1. INTRODUCTION 

The empirical relationship between masses of supermas¬ 
sive black holes (BHs) and the stellar velocity disper- 
sion or integrated stellar luminosity of their host galaxies 
(lKormen dvll l993HMagorrian et al.lll998HEe rrarese & Merritt! 

I2000t iGebhardt et al.l l2000t iGiiltekin et al.l 120091) has been 
taken as an indication of an evolutionary connection between 
the black hole and host galaxy. The discovery of this relation¬ 
ship has motivated the development of a new generation of 
merger-driven galaxy evolution models in which coordinated 
black hole accretion and star formation in the host galaxy 
lead to the observed relation between thei r masses (e.g., 
\^lonteri_etalj_|2003 : jpj_Matteo et ^ 120051 iHopkins et al.l 
2006t I Volonteri et al. I2015L and references thereinl. The the¬ 


oretical models are bolstered by observations of dual ac¬ 
tive galactic nuclei (AGNs), with kiloparsec-scale separa¬ 
tions, in interacting and merging galaxies. These can be 


detected directly via imaging (e.g.. 

Junkkarinen et al.l 1200 It 

Komossa et al.l I2003t 

Ballo et al.l 

2004 

t Comerford et al.l 

2009at lEabbiano et al.l 

20111 IKossetal. 

120111 1201211 and 


indirectly via their double-peaked narrow emission lines 

(IComerford et al.ll2009bt IWang et ^l2009t ISmith et ^1201 Ot 
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The formation of supermassive black holes binaries (SB¬ 
HBs) during the late stages of the mergers is one outcome 
of hierarchical galaxy evolution. In the e volutionary sce¬ 
nario of SBHBs, first described in detail by iBegelman et alJ 
(Il980h . the orbital decay of close, bound binaries slows 
down dramatically or stalls at separations of order 1 pc 
because of “loss cone” depletion (the loss cone is the re¬ 
gion in stellar phase space containing stars that can scat¬ 
ter off the SBHB and drain its orbital angular momen¬ 
tum). This slow-down has been termed the “last parsec prob¬ 
lem.” However, modem calculations have shown that the 
SBHB orbital decay may proceed unimpeded becaus e the loss 


20051 iMilosavlievic & Merritt! 

2003 

), the influence of non- 

spherical potentials (e.e.. lYul 

2002 

; iMerritt & PoonI 120041 

Berczik et all 12006 

; iKhanetalJ 120131 IVasiliev & MerrittI 


20131 IVasiliev et al.1120147 or interactions between the bi 


narv and a reservoir of gas ("e.g.-lArmitage & Nataraianll2002 

Escala et al.l 

2004lDotti et al.120071 2009llCuadra et al.ll2009 

Lodato et al. 

l2009llParriseta1.||2Q15a|i. 




Observational evidence for close, bound SBHBs has 
proven quite elusive. The two noteworthy candidates are 
CSO 0402-1-379, with a projected separat ion of approximate! 

7 pc (imaged by radio int erferometry iManess et alJ 1200- 
iRodriguez et alJl2006L l2009h and OJ 287, where outbursts in 
the light curve occurring on a ^12 year period h ave been 
explained in this context (e.g., IValtonen et alJUOlX and ref¬ 
erences therein). The identification of such systems is par¬ 
ticularly difficult because these small separations are not re¬ 
solved at cosmological distances. Yet, finding SBHBs is im¬ 
portant because it would validate the galaxy evolution sce¬ 
narios mentioned above, because SBHBs have been invoked 
to explain a variety of other astro nomical observations (see 
summary in lEracleous et al.l|2012L hereafter Paper I), and be¬ 
cause they are the progenitors of low-frequency gravitational 
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wave sources potentially detec tabl e by Pulsar Ti ming Arrays 
(e.g.. lArzoumanian et al.llMTil see lSesan al|2015|f or a review) 
and f uture space-based instruments (see lAmaro-Seoane et aU 
120121 for a review). 

Recent observational searc hes for sub-pa rs ec SB HBs, 
motivated by the wor k of iKomossa et al.l (l2008h and 
iBoroson & Lauej (l2010h . have turned up many new candi- 
dates. An in-dep th review of the field can be found in 
iBogdanovid (|20 1 5h and a discussion of an array of techniques 
to identify SBHBs , primarily via their emission-lines, is given 
in iPopovid (l2012h . The widespread approach has been to 
make a n analogy with single-line spectroscopic binary stars 
(see e.g. lGaskelll 19961 for a review of this concept) and search 
for indications of orbital motion in the periodic velocity shifts 
of the broad emission lines relative to the narrow lines. In this 
spirit, a number of groups have searched for quasars from the 
Sloan Digital Sky Survey (SDSS) whose broad H/3 or Mg II 
emission lines are displaced from their systemic redshifts by 
several hundred to a few thousand km s“* and/or show sub¬ 
stantial radial v elocity variations about their systemic red¬ 


shifts (Paper I: iTsalmantza et al.l 1201 It 

Decarli et al.l 12013t 

iShen et alJl2013tlLiu et al.ll2014ll.Iu et all 

2013h. In the work- 


ing hypothesis adopted for these searches, the SBHB sits in a 
circumbinary disk and the secondary, which can more easily 
access gas on the inner edge of this disk (see, for example, 
lArtymowicz & Lubowlll99a iHayasaki et alJl200^ . accretes 
at a much higher rate than the primary and alone has a broad¬ 
line region (BLR). In the picture emerging from the observa¬ 
tions, SBHB candidates that have been extensively tested and 
are still widely held to be viable are rare, but initial studies 
of pairs of spectra, taken years apart, suggest that in some 
cases observed velocity shifts of the broad emission lines are 
consistent with orbital motion (e.g.. Paper I lShen et'^l2013l : 
Liu et aljr2013l : lLiu et al.ll2014^ . It is important to note that the 
interpretation of these measurements is complicated by many 
significant caveats, as discussed in Paper I, because the broad¬ 
line shifts are not unique signatures of orbital motion. 

Obtaining more robust constraints on the nature of the 
SBHB candidates via the radial-velocity technique can be ac¬ 
complished by long-term spectroscopic monitoring. If the or¬ 
bital periods are of ord e r a century, as noted by a numb er of 
authors (e.g.. lYull2002L lLoebll201(]l : iDecarli et al.ll2()T^ Pa¬ 
per I), we may be able to detect small but systematic radial ye- 
locity yariations oyer the course of 1-2 decades. If the orbital 
periods happen to be of order a few decades, we may be able 
to monitor a complete orbital cycle in our lifetimes. Thus, we 
haye been carrying out a spectroscopic monitoring program of 
the sample from Paper I and refining our methodology for an¬ 
alyzing the data. In this paper, our primary goal is to present 
the spectra we haye obtained since Paper I and a refined anal¬ 
ysis aimed at (a) yerifying some of the results presented in 
Paper I and (b) quantifying the uncertainties in our measure¬ 
ments of spectroscopic properties. We also take adyantage 
of this opportunity to pursue a secondary goal: explore the 
yariability of the continuum and broad H/3 line fluxes of our 
SBHB candidates and compare them to those of a sample of 
typical quasars of similar luminosity and redshift. In a forth¬ 
coming paper we concentrate on the study of long-term radial 
yelocity yariations, where we will adopt lessons learned here 
about broad line yariability. 

The continuum emission in quasars is known to yary, the 
amplitude of fluctuations is typically on the order of 10-30%, 
increasing towards longer time scales and eyentually leyel- 
ing off at periods longer than a few years (as inferred from 


long; 


1993 


t erm light curyes and structure functions: e.g.jHawkins 


iGiyeon et al.l 119991 : ISesar et akl l2006t iMacLeod et al 


2012L and references therein). Howeyer, the mechanism pro¬ 


ducing the offset emission lines in the absence of a SBHB 
may leaye its imprint on the flux yariability as well. Recoiling 
black holes, where the supermassiye black hole and BLR are 
kicked out of the center of the host galaxy following a merger, 
can also produce offset broad lines and would presumably 
display yariability of the continuum and broad emission line 
flux that is consistent with a typical quasar. Although simula¬ 
tions suggest that the accretion history of recoiling black holes 
will differ significantly from ordinary AGN oyer the course of 
their 10-100 Myr lifetimes (iBlecha et al.ll201 it iGuedes et alJ 
ISoOb, on time scales of decades or less the yariability proper¬ 
ties should be normal. It is not completely clear what might be 
expected of flux yariability in a SBHB. Simulations of mas- 
siye black hole binaries wi th gaseous disks generally And yari- 
able accretion rates (e.g . ,jHagasaki_etaLlj2b07yCuadraeg^ 
2(M iNoble et all 12012: iRoedig et alJl2012L iShi et alJl2012L 


Bode et al.l2012LlD’Orazio et alJ201^lFarris et al.l201^ . and 

calculations suggest t hat the luminosity a nd accretion rate can 
remain proportional (iFarris et al.ll2015bh . but beyond this the 
results diyerge. The time scale for yari ability can range from 
0.5 to 10 times th e binary period (e.g.. lD’Orazio et ffll2013L 
IFarris et ani2014f) . Obseryational searches haye taken adyan¬ 
tage of the possibility that continuum yariability may be syn¬ 
chronized with orbital motion of the binary and some SBHB 
candidates haye been claimed based on their light curye (e.g., 
iGraham et al.ll2015t iLiu et alJl2015h . In these cases the or¬ 
bital periods are yery short, amounting to a few years or less. 
Taken at face yalue, such systems would represent the last 
stage of the eyolution of the SBHB, just before the merger 
(in comparison, the candidates found yia offset broad emis¬ 
sion lines haye orbital periods of order a few hundred years 
and, if their orbits were to decay by grayitational radiation 
alone, tit would take them ^ 10'° years to reach periods of 
a few years). Howeyer, the ubiquity of red noise in quasar 
power s pectra introduces subs tantial uncertainty into this ap¬ 
proach dMacLeod et alJiTolOh . With these considerations in 
mind, studying the flux yariability properties of objects with 
offset broad emission lines can be yaluable and informatiye, 
whether the line offsets are indicatiye of binaries, recoiling 
black holes, or an extreme population of indiyidual quasars. 

In Section |2] we summarize the sample selection and de¬ 
scribe the acquisition of new data and the construction of a 
comparison sample of typical quasars. In SectionOwe assess 
the ensemble yariability of the binary candidates and com¬ 
parison sample, both by considering the changes in the in¬ 
tegrated H/3 and optical continuum fluxes and by comparing 
structure functions for the two samples. We discuss our re¬ 
sults in the context of preyious work and the binary hypoth¬ 
esis in Section IH and also consider the constraints implied 
by the obseryed yariability on the properties of the binaries. 
We summarize this inyestigation in Section |5] Additionally, 
Appendix lAl proyides spectral measurements that haye been 
made for the sample of SBHB candidates but are not utilized 
in this work. Throughout this work, we adopt a cosmology of 
Ho = 73 km s“* Mpc“', Ha = 0.73, and Q,„ = 0.27. 





















































































Hux Variability of Supermassive Black Hole Binary Candidates 


3 



Figure 1. Distribution of redshifts for the SBHB candidates and the com¬ 
parison sample of typical quasars from the SDSS. The two distributions are 
clearly not identical but they span approximately the same range and have 
approximately the same median. 

2. DATA 

2.1. The Sample of Supermassive Black Hole Binary 
Candidates and New Spectra 

The selection process and spectroscopic properties of the 
sample of SBHB candidates are described in detail in Pa¬ 
per I. In summary, the sample was selected via spectral prin- 
ci pal component analysis ( following the technique described 
in iBoroson & Lauerl 120 id) to identify a population of ob¬ 
jects from the SDSS whose broad H/3 profiles do not re¬ 
semble those of typical quasars. We then vetted the initial, 
automatically-selected sample by visual inspection, discard¬ 
ing, for example, clearly double-peaked profiles (which are 
likely produced by an accretion disk around a single super¬ 
massive BH, e.g., |Eracleo us et al.lll997HEracleous & HalpernI 
120031 lEracleous et al.ll2009h . The final sample contains 88 
objects having a median redshift of 0.32, with a full range 
of 0.077 < z < 0.713, as shown in the top panel of Figure [T] 
Absolute V-band magnitudes, after Galactic extinction cor¬ 
rections but not K corrections, range from -21.25 to -26.24, 
with a median value for the sample of -23.0^ 

By virtue of the sample selection, the initial spectroscopic 
observations come from the SDSS. Paper I presents additional 
spectra for many objects in the sample from an array of tele¬ 
scopes and instruments. For this work, we use all the spectra 
from Paper I, plus new observations that are presented here 
for the first time. In Table [1] we give a summary of the new 
spectroscopic observations. For each of the 212 new spectra 
we give the observation date, telescope and instrument used, 
exposure time, signal-to-noise ratio (S/N) in the continuum 
near the line of interest, and rest-frame wavelength coverage. 
For completeness, we also include the redshift of each object, 
the apparent and absolute V magnitude, and the Galactic V- 
band extinction (see the footnotes to Table [T] and section 4. 1 
of Paper I for further details on the target properties) and we 
also list objects for which no new observations are presented. 

* We use “quasar” and AGN interchangeably in this work, even though 
some objects in this sample do not meet the strict absolute magnitude crite¬ 
rion for quasars. 
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Figure 2. For every unique pair of observations for the SBHB candidates, 
the distribution of rest-frame time intervals (top). Am = -2.51og(/i//2) mea¬ 
sured from continuum flux density ratios at 5100 A (bottom right), and Am 
at 5100 A as a function of rest-frame time interval between observations (bot¬ 
tom left). The Am distribution has a median of—0.05 and a standard deviation 
of 0.43. 



Figure 3. For every unique pair of observations for the SBHB candidates, 
the distribution of rest-frame time intervals (top). Am = —2.51og(/i 7/2) mea¬ 
sured from integrated H/3 flux ratios (bottom right), and the integrated H/3 
Am as a function of rest-frame time interval between observations (bottom 
left). The Am distribution has a median of—0.03 and a standard deviation of 
0.59. 


For reference, we note that in Paper I we presented the origi- 
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Table 1 

Log of Spectroscopic Observations [Abridged] 


Object Name 

SDSS J 
(1) 

(2) 

mv*’ 

(mag) 

(3) 

Av*’ 

(mag) 

(4) 

(mag) 

(5) 

No. 

obs.^ 

(6) 

Observation 
Date (UT) 

(V) 

Instr. 

Code' 

(8) 

Exposure 

Time 

(s) 

(9) 

S/m 

(10) 

Rest-Frame 
Wavelength 
Range (A) 
(11) 

001224.03-102226.2 

0.2287 

17.07 

0.127 

-23.27 

5 

2010/10/12 

H2 

1798 

27 

3481-5911 







2014/08/28 

Ab 

1560 

41 

2848^476 







2014/08/28 

Ar 

1560 

27 

4197-8013 

002444.11+003221.3 

0.4024 

16.85 

0.083 

-24.86 

4 

2014/08/29 

Ab 

1260 

49 

2495-3921 







2014/08/29 

Ar 

1260 

52 

3691-7033 

015530.02-085704.0 

0.1648 

16.84 

0.080 

-22.65 

5 

2011/08/31'* 

H2 

2587 

10 

3681-6238 







2014/08/28 

Ab 

1260 

28 

3004-+1721 







2014/08/28 

Ar 

1260 

10 

4426-8451 


^ The redshift of the object as measured in this paper from the peak wavelength of [O III] A5007 line. 

^ The apparent V magnitude, determined from the SDSS PSF magnitudes, as described in §4.1 of Paper I. 

The Galactic visual extinction, taken from lSchlegel, Finkbeiner, & DavijU998l) 

^ The absolute V magnitude, computed as described in jOof the text. 

® The total number of observations, including those reported in Paper I. 

^ The telescope and instrument configuration, as described in Table 

® The signal-to-noise ratio (S/N) in the continuum near the line of interest. When the spectmm includes the H/9 line we give the S/N in the continuum 
near this line, at 4600 A. If the spectnam includes only the Mg II A2800 line, we give the S/N in the continuum near this line, at 2900 A. Some of the 
APO spectra cover the continuum between the above lines but neither of the lines themselves; for these spectra we report the S/N in the continuum at 
3600 A. 

^ Average of two spectra taken within a few days of each other. The exposure time is the sum of the individual exposure times and the date is the 
“median” of the dates of the two observations. 


Table 2 

List of Telescopes and Instruments 


Instrument 

Code 

Observatory, Telescope, and Spectrograph 

Spectral Elements 

Resolution^ 

(A) 

MM 

MDM, Hiltner 2.4m, MODSPEC spectrograph 

600 mm“* grating, K'Oslit 

3.4 

MO 

MDM, Hiltner 2.4m, OSMOS spectrograph 

VPH grism (704 mm^'), 172 slit 

3.8 

K 

KPNO, Mayall 4m, Ritchie-Cretien spectrograph 

BL420 grating (600 mm“*), 175 slit 

2.4 

Ab 

APO, ARC 3.5m, Double Imaging Spectrograph 

blue arm: 400 mm“^ grating, 1''5 slit 

6.3 

Ar 


red arm: 300 mm“^ grating, 1''5 slit 

7.3 

H2 

Hobby-Eberly Telescope, Low-Resolution Spectrograph 

G2 grism (600 mm“*), l75 slit 

5.6 

H3 


G3 grism (600 mm“^), VJ5 slit 

4.9 


^ The spectral resolution (FWHM) at 6400 A (the median wavelength of H/9 for the redshift distribution of our targets). 


nal SDSS spectra for our 88 targets plus 108 followup spectra. 
We also tabulate the total number of observations per object; 
with the addition of the new spectra, all objects in the sample 
have at least two observations and most objects have at least 
three observations. Time intervals between spectra in the ob¬ 
served frame range from weeks to over eleven years. The 
distribution of rest-frame time intervals is shown graphically 
in Figures|2]and[3 

The new spectra presented here were obtained with the Hilt- 
ner 2.4m telescope at the Michigan-Dartmouth-MIT obser¬ 
vatory (MDM), the Mayall 4m telescope at Kitt Peak Na¬ 
tional Observatory (KPNO), the Astrophysical Research Con¬ 
sortium (ARC) 3.5m telescope at Apache Point Observatory 
(APO), and the 9.2m Hobby-Eberly Telescope (HET). The 
instrument configurations used to obtain the spectra are listed 
in Table |2] along with the spectral resolution attained by each 
instrument configuration. The spectra were reduced and cal¬ 
ibrated as described in section 5.1 of Paper I. Since accu¬ 


rate wavelength calibration is important for our purposes, we 
summarize the wavelength calibration process here. We used 
20-60 arc emission lines to derive the wavelength solution 
for each observing run. This solution consists of a polyno¬ 
mial of order five or less, connecting the detector pixel num¬ 
ber to the wavelength. We required that the standard devi¬ 
ation of fit residuals was 0.1 pixel or less. Thus, the rela¬ 
tive velocity scale of each spectrum at 6400 A (the median 
wavelength of for our quasars) is good to 5 km s“^ for 
MDM spectra, 7 km s“^ for KPNO spectra, 11 km s“' for 
APO spectra, and 9 km s“' for HET spectra. We used the 
[O l] A5577 line in the night-sky spectrum, recorded simul¬ 
taneously with and on the same detector as the quasar spec¬ 
trum, to set the absolute wavelength scale of our new spec¬ 
tra. After applying heliocentric velocity corrections and con¬ 
verting the wavelength scale to vacuum values, we aligned 
the quasar’s [O III] AA4959, 5007 doublet in each new spec¬ 
trum to that of the original SDSS spectrum using the cross- 
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Figure 4. The distribution of continuum luminosity for the SBHB candi¬ 
dates and comparison sample. Although the distributions ai'e not identical, 
the samples cover a similar, relatively narrow range in luminosity which al¬ 
lows a reasonable comparison independent of luminosity effects. 


correlation method described in Paper I. The final uncertainty 
in this alignment translates to an uncertainty in relative veloc¬ 
ities between the broad H/? lines observed at different epochs. 
As noted in Paper I, this uncertainty is smaller than 25 km s“* 
for all objects and smaller than 7 km s“' in 80% of objects in 
our sample. 



Figure 5. For every unique pair of observations for the comparison sample, 
the distribution of rest-frame time intervals (top), Am = -2.51og(/i//2) mea¬ 
sured from continuum flux density ratios at 5100 A (bottom right), and Am 
at 5100 A as a function of rest-frame time interval between observations (bot¬ 
tom left). The Am distribution has a median of 0.03 and a standard deviation 
of 0.22. 


2.2. A Comparison Sample of Typical SDSS Quasars 

In order to characterize the variability properties of typical 
quasars, we generated a comparison sample of SDSS quasars 
with multiple spectroscopic ob servations in the seve nth data 
release (DR7) quasar catalog (ISchneider et al.ll2007^ . Typi¬ 
cally, each object is observed two to three times, although sev¬ 
eral objects have as many as ten observations. Starting with a 
list of objects with multiple observations, we selected those in 
the redshift range 0.08 < z < 0.7, which matches that of our 
SBHB candidates. To ensure our ability to make reliable mea¬ 
surements from the spectra, we also requir ed that S/N >20 i n 
the H/3 region (using the S/N values from IShen et al.ll^l Ibh . 
These criteria resulted in a comparison sample of 212 quasars 
that we analyze in an identical way to the sample of SBHB 
candidates throughout this work to facilitate a fair compari¬ 
son free of any methodological biases. 

All observations of quasars in the comparison sample were 
carried out with the SDSS spectrograph through fibers of di¬ 
ameter 3". The resulting spectral resolution was 2.7 A, which 
is slightly better than that achieved in followup observations 
of the sample of SBHB candidates (see Table|2]i. Differences 
in spectral resolution between SDSS spectra and followup 
spectra is inconsequential for our purposes since all the lines 
of interest, including the narrow H/3 and [O III] lines, are al¬ 
ways broader than the instrumental resolution. 

Given that quasar variability properties are known to de¬ 
pend on luminosity and redshift, it is important that the red- 
shifts and luminosities (primarily the luminosities) of the 
comparison sample span the same range as those of the SBHB 
candidates. The redshift distribution of the quasars in the 
comparison sample is shown in the bottom panel of Figure [T] 
The redshifts span the range 0.08 < z < 0.68, with a me- 



Figure 6. For every unique pair of observations for the comparison sample, 
the distribution of rest-frame time intervals (top). Am = —2.51og(/i 7/2) mea¬ 
sured from integrated H/3 flux ratios (bottom right), and the integrated H/3 
Am as a function of rest-frame time interval between observations (bottom 
left). The Am distribution has a median of 0.02 and a standard deviation of 
0.28. 


dian value of 0.38, very similar to the binary candidates de- 
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scribed in Section ITTl The monochromatic luminosity dis¬ 
tributions of the binary candidates and t he c omparison sam¬ 
ple are shown in Figure |4] (see Section 12.31 below for a de¬ 
scription of the luminosity measurements). The binary can¬ 
didates span the range 43.27 < log [ALa( 5100A)/erg s“'] < 
44.68 with a median value of 43.75 and a standard deviation 
of 0.31. The comparison sample spans the range 43.44 < 
log [ALa( 5100 A)/erg s“*] < 45.17 with a median value of 
43.88 and a standard deviation of 0.33. The shapes of the red- 
shift and luminosity distributions of the SBHB candidates and 
comparison quasars are not identical, of course, as indicated 
by the Kolmolgorov-Smirnov (KS) test, which gives probabil¬ 
ities that the redshift and luminosity distributions were drawn 
from the same parent population of 0.002 and 0.006, respec¬ 
tively. 

The distribution of rest-frame time intervals is shown in 
the top panels of both Figure |5] and Figure |6l The me¬ 
dian, observed-frame time interval between observations of 
the same object in this sample is approximately 225 days, 
with concentrations of time interval between observations at 
less than 100 days and approximately 300 days (these are not 
obvious in relevant histograms because a very fine binning 
would be required to resolve them, but the concentrations are 
evident in the associated scatter plots). The minimum time 
interval in the observed frame is less than a day, and the max¬ 
imum is nearly eight years. 

2.3. Spectral Decomposition and Emission-Line 
Measurements 

In order to isolate the broad H/3 line and quasar continuum, 
we perform a spectral decomposition that deblends the quasar 
continuum, optical Fell, H/3, and [Olll] emission compo¬ 
nents. All line and continuum f itting is performed using the 
IRAF task SPECFIT (iKrissll 19941) in three steps. This process 
is qualitatively similar to the approach of Paper I and was 
adopted because, in objects where the line profiles are partic¬ 
ularly complex, it is more effective at characterizing the red 
wing of H/3 under the [Olll] emission t han methods that f it 
all the components simultaneously (e.g.. IShang et al.l 120051) . 
In fact, attributing various components of the fit to different 
physical regions is not advisable, but it is unavoidable in this 
case, and we feel that this method minimizes the danger of 
misrepresenting the data. 

In the first step, we fit and subtract the optical continuum, 
which includes contributions that can be characterized by the 
combination of a featureless powe r law of the form f \ oc A~° 
and the optical Fell template of iVeron-Cettv et'ffl (12004ft . 
We adopt two suites of two initial guesses, where the initial 
guesses differ in the properties of the starting power law. The 
first suite includes only the power-law continuum and Fe II 
template. The second suite adds a Gaussian to account for 
He II A4686 emission and was only selected when the He II 
emission was visually evident in the spectrum. The fit is opti¬ 
mized in three windows at 4150^200,4435^700, and 5100- 
5700 A and the best fit is selected by applying the Ba yesian 
Information Criterion (hereafter BIC, ISchwartzllI978ft . Al¬ 
though it is appropriate to minimize the statistic when per¬ 
forming the fit, the reduced is ineffective at choosing be¬ 
tween models because a fit with more free parameters will 
always perform better at this test. The BIC provides an eval¬ 
uation of the goodness of fit, similar to that provided by the 
reduced x^, but penalizes the model based on the number of 
free parameters, and thus evaluates whether the addition of 
new free parameters truly improves the fit. In a minority of 


SDSSJ133432.35+171147.0 



Rest Wavelength [A] 


Figure 7. An example of the spectral decomposition for SDSS 
J133432.35+171147.0. The red line shows the total model, including the 
power-law continuum, Fe II template (which is weak in this object), the nar¬ 
row [O III] lines, and the broad and nan'ow H/3 lines. The individual line 
profiles are shown in blue. 

cases, the optical continuum was not well fitted by the above 
procedure; For example in some cases where data were miss¬ 
ing in one of the fit windows. In these cases, manual adjust¬ 
ments were made to the initial guess or wavelength range of 
the fit. Finally, continuum-subtracted spectra were generated 
by subtracting the best-fitting optical continuum model from 
the data. 

In the second step, we characterize the [O III] AA4959, 5007 
emission lines and subtract them. We parameterize the profile 
of narrow emission lines in the spectrum based on the profile 
of the [O III] A5007 line. In order to isolate that line, we fit 
a low-order polynomial to the local continuum (i.e. the red 
wing of the broad H^ line) in manually selected wavelength 
windows. The polynomial is temporarily subtracted and the 
remaining [Olll] A5007 line is fitted with two Gaussians to 
adequately account for any possible profile asymmetry (no 
physical meaning is attributed to the individual Gaussians). 
The [O III] A4959 line profile is then taken to be identical to 
that of [O III] A5007, but appropriately shifted in wavelength 
and scaled down in flux by a factor of 1/3. With the [Olll] 
profiles thus characterized, the low-order polynomial is added 
back to the data and the [O III] doublet is subtracted to isolate 
the H^ profile. 

In the third and final step, we decompose the H/3 profile. In 
typical quasars, the H/3 profile can be adequately represented 
by a combination of four Gaussian components. Two of the 
Gaussians account for the narrow H/3 line and are constrained 
to create a profile identical in shape and tied in wavelength to 
the [Olll] A5007 line, although the flux of the profile is al¬ 
lowed to vary and even to be zero if this provides the best fit. 
The other two Gaussians are allowed to vary freely to charac¬ 
terize the broad H/3 line. In cases where these profiles are very 
complex, an additional Gaussian is necessary to adequately 
characterize them. These cases are identified by visual inspec¬ 
tion, and are outliers in the distribution of reduced x^ when 
we fit them with only two Gaussians (typically, jv > 16). 
Including the third Gaussian results in a decrease in the BIC 
of > 10, indicating that the additional free parameters associ¬ 
ated with the third Gaussian significantly improve the fit. An 
example of the fitting results is shown in Figure|7] 

We measure spectral properties from the parametric models 
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obtained for the H/3 and [Olll] A5007 lines in the SDSS spec¬ 
tra and give the results in Tables[3and|4] The total line fluxes 
are obtained by integrating the emission-line model where the 
model is above 1% of the peak flux density of the line. We 
verify visually that this scheme includes the entire line flux 
(see Ar)r)endix lA.2l for a detailed discussion). The continuum 
luminosity is measured at 5100 A from the power-law com¬ 
ponent of the continuum fit and is included in Table 01 The 
continuum luminosity distributions for the SBHB candidates 
and comparison sample are shown in Figured The luminos¬ 
ity range in both samples is relatively narrow, spanning only 
about an order of magnitude. 

Using the spectral decompositions, we also determine a 
suite of parameters to characterize the H/3 and [Olll] A5007 
line profiles. Although not all of these quantities are needed 
for the variability analysis carried out in this paper, we in¬ 
clude them in Tables [3 and 0] for future reference. Following 
Paper I, the first four central moments of the line profiles, p„, 
are calculated from fi„ = where A; and fi 

i 

are the discrete wavelength and flux densities, the normaliza¬ 
tion constant K is given by 1 /K = and (A) = KY^ifi 

i i 

is the first moment (i.e. the centroid) of the line profile. 
From these measurements we can derive the skewness coeffi- 
cient, i , which describes the degree of asymmetry 

of the profile. Symmetric profiles have i = 0 and positive val¬ 
ues of the skewness coefficient indicate that the line leans to 
bluer wavelengths (i.e., the red wing is more extended than 
the blue). The kurtosis coefficient, k = ^* 4 /^ 2 ’ describes the 
shape of the line profile, with boxy profiles having smaller 
values of k. We also calculate the Pearson skewness coef- 
ficient, p = {{\) - \,„)/pY. In this case, A,„ is the median 
wavelength corresponding to the point where the area of the 
line profile is split in half. Note that this definition matches 
that given in Paper I, but the values given in Table 4 of Paper I 
had the sign reversed, inadvertently. We correct this mistake 
in this paper. Positive values of the Pearson skewness coeffi¬ 
cient also describe blue-leaning line profiles. 

An accurate accounting of the uncertainties associated with 
the spectral properties must include both the effect of noise in 
the spectrum as well as uncertainties incurred during the in¬ 
teractive fitting procedure. To estimate these uncertainties we 
carry out Monte Carlo simulations of the spectral decomposi¬ 
tion. First, we select five spectra that span the range of S/N 
present in the sample of SBHB candidates, making an effort to 
include instances of qualitatively different spectral properties. 
We calculate the S/N within ±25 A of the peak of the broad 
H/3 profile after continuum subtraction. For each of these five 
objects, we generate fifty realizations of the best-fitting model 
spectrum by adding Gaussian noise with a constant amplitude 
equal to that measured at the peak of the broad H/3 line. Each 
synthetic spectrum is then run through the spectral fitting pro¬ 
cedure (including the interactive continuum placement and 
line fitting) and all of the spectral properties are measured 
by fitting the synthetic spectrum in the same manner as the 
observed spectrum (see description above). We adopt the re¬ 
sulting standard deviation in the distribution of each measured 
property as the uncertainty for that measurement at that S/N. 
As an example of the outcome of this experiment, we show 
in Figure 0 the distribution of H/3 FWHM values for differ¬ 
ent values of the S/N. Next, we fit a simple function of the 
form a =A(S/N)~°‘+B to the relation between the S/N and 
corresponding uncertainty, as seen in the example of Figure|9l 



-1000 -500 0 500 1000 
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Figure 8. Distributions of the FWHM of H/3 measured from the 50 synthetic 
spectra generated for 5 different objects. Each panel represents a different 
object and is labeled with the value of S/N measured from the observed 
spectrum. The standard deviations of these distributions yield the FWHM 
uncertainties as a function of S/N. 



HP S/N 

Figure 9. Uncertainties in the FWHM taken from the widths of the distribu¬ 
tions in the previous figure versus S/N. These are fitted with a function of 
the form a = A{S/Ny^ +fi, where A = 2341.33, B = -103.40, and a = 0.71 
in this case. 

We find that the uncertainty usually scales as either (S/NY 
or (S/NY^^. In a handful of cases, we fit the fractional uncer¬ 
tainty (relative to the mean of the distribution of realizations) 
rather than the absolute uncertainty as a function of S/N be¬ 
cause we find that this relationship is much more clear. Using 
these relationships we estimate uncertainties on the measured 
quantities for the objects in the sample. For the FWHM and 
peak wavelength shift of [Olll], we find a negligible spread 
in the values measured from the synthetic spectra, likely be¬ 
cause the line profiles have strong peaks and few data points, 
therefore we set the following upper limits to these uncertain- 
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ties: 70 km s“' and 30 km s“*, respectively. In Appendix lAl 
we investigate the sources of uncertainty in detail, including 
the uncertainty arising from the choice of method for decom¬ 
posing the spectra. The results presented in the tables of this 
paper do not include the uncertainty from the choice of de¬ 
composition method. 

The primary objective of the multi-epoch spectroscopy of 
the SBHB candidates is to search for radial velocity varia¬ 
tions, which will be the topic of an upcoming paper. As a re¬ 
sult, while particular care was taken with the wavelength cali¬ 
bration, the absolute flux calibration can be uncertain by up to 
a factor of two, largely because of variable seeing. To bypass 
this uncertainty, we assume that the integrated [O III] lumi¬ 
nosity does not change over the time scales spanned by our 
observations. To support this assumption, we verify that the 
narrow-line regions (NLRs) of our targets have a smaller an¬ 
gular size than the spectrograph aperture (circular fibers with 
a radius of 1"5 in the case of SDSS spectra and a slit with 
a width of 1 "5 in the case of follow-u p spectra). To do this , 
we use luminosity-based prescription of iBennert et al.l (|2002), 
which gives a measure of the outer extent of the NLR (Rnlr, 
defined as the radius of the ring aperture where the [O III] flux 
in their images drops to 3a above the background and typi¬ 
cally enclosing 98% of the detectable emission) as a function 
of the [Olll] luminosity. The scatter in this relationship is 
0.14 dex in \og{RffLR), so in order to obtain an upper limit on 
the size of the NLR in our objects, we add this to the val¬ 
ues calculated from the prescription. For maximum sizes of 
the NLRs in our sample, we use the [O III] luminosity from 
the SDSS spectrum for each object (because these have the 
best flux calibration) and obtain 800pc < RNLR,mca < 14600pc 
with a median value of 2700 pc and corresponding angular 
sizes of 0."12 < a < 1 ."46 with a median value of 0."42. Based 
on this calculation, none of the objects in the sample have 
maximum NLR sizes larger than L"5. The angular sizes that 
we estimate are a lso consistent with expecta tions from more 
recent wo rk (e.g..lHainline et al.ll2 013Ll201 4^. G iven that the 
images of iBennert et akl (^2002^ and lYoung et al.l (1201 4ft show 
that the NLR emission is actually concentrated far inside of 
the outermost extent that is measured by Rnlr, we determine 
that there should be no variability as a result of a small size 
or orientation of the spectroscopic slit. Intrinsic narrow-line 
variability in response to changes in the ionizing continuum 
can take place on time scales of years (e.g., NGC 5548 which 
demonstrates about 20 % variation over 20 ye ars, suggesting 
that t he NLR extends down to a few parsecs, iPeterson et al.l 
1201 3ft . Because the NLR is so large in the binary candi¬ 
dates and comparison sample, which are more luminous than 
NGC 5548, we do not expect intrinsic variability on the time 
scales of our observations. Thus, whenever we compare two 
spectra from this point forward, we scale all flux densities and 
integrated line fluxes so that the integrated [O III] flux of the 
later observation matches that of the earlier one. 

3. ANALYSIS 

3.1. The Distribution of Magnitude Differences 

In this work, we consider the rest-frame, ensemble flux 
variability in the continuum and H/3 emission line. Because 
each object in the SBHB candidate sample has only been ob¬ 
served a handful of times, the light curves are poorly sampled 
and cannot provide a useful characterization of the variability 
of individual objects. Therefore, we consider the variability 
of the sample of SBHB candidates as a whole by means of the 


structure function. In this approach, we consider every possi¬ 
ble pair of spectra for each object to obtain a measure of the 
dispersion in the magnitudes, where Am = -2.5 log(/i 7 / 2 ), at 
different time scales for the entire sample. 

For the SBHB candidates, the right panels of Figures |2] and 
[3show the distributions of continuum and integrated H/3 mag¬ 
nitude differences, respectively, of unique pairs of observa¬ 
tions. Although the distribution is centered near zero there 
are noticeable wings that are more pronounced for H/3 than 
the continuum. The distribution of continuum magnitude dif¬ 
ferences has a median of-0.05 and standard deviation of 0.43 
while the distribution of H/3 magnitude differences has a me¬ 
dian of-0.03 and a standard deviation of 0.59. 

The distribution of magnitudes does not take into account 
the time interval between observations, which is known to 
play a role in dete rmining the amplitude of flux variability 
in typ ical quasars (IVanden Berk et all 120041 : iMacLeod et alJ 
l2012ft . Figures |2] and [3 also show the continuum and inte¬ 
grated H/3 magnitude differences, respectiyely, for the SBHB 
candidates as a function of the rest-frame time delay between 
the obseryations. It is readily apparent that a substantial frac¬ 
tion of obseryations are separated by 1000-3000 days. This 
is a direct result of the obserying strategy we adopted, where 
the goal was to eyaluate the binary hypothesis based on the 
H^ radial yelocity yariability. It is also noteworthy that the 
distributions of continuum and integrated H/3 magnitude dif¬ 
ferences are fairly similar. 

3.2. Ensemble Variability 

In order to quantify the ensemble yariability of the SBHB 
candidates as a function of the time interyal between obser¬ 
yations, we compute the structure function, which proyides 
a measure of the ayerage relatiye change in the flux of the 
objects in the sample. There are multiple definitions in the 
literature, we adopt two here. The first yersion and the cor¬ 
responding uncertainty, from IVanden Berk et al] (12004ft . are 
defined as follows: 

SF{At) = ^'!^{\Am\Y- {a\^f) (1) 

^ ^^■ ( 2 ) 

We eyaluate the magnitude difference Am = 
-2.5 log(/i// 2 ) between two epochs (rest-frame time 

interyal At), after scaling the later spectrum to the same 
[Olll] flux as as the former spectrum. The quantity (TAm 
represents the formal measurement uncertainty in Am and 
arms.Am and CTraij o-A,,, are the standard deyiations of the 
Am and ctahi distributions. The structure function is then 
calculated in time interyal bins. 

The sec ond yersion of the struc ture function we consider is 
defined by IMacLeod et al.1 (12012ft as 

SF,Qu(At) = 7[0.74/eR(Am)]2-aL„ (3) 

where IQR is the interquartile range (full width uLquarter 
maximum) of the cumulatiye Am distribution and aj^^ is the 
median of the squared Am measurement uncertainties. The 
uncertainty in this structure function is simulated by bootstrap 
resampling of the distribution of Am 1000 times, calculating 
the corresponding yalue of the interquartile range for each it¬ 
eration, and finally eyaluating the width of the resulting dis¬ 
tribution of SFiqr. 
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Table 3 

Spectral measurements for broad H/3 [Abridged] 


Object 

Name 

SDSSJ 

(1) 

Obsev. 

Date 

(UT) 

(2) 

FWHM 

(km/s) 

(3) 

FWQM 

(km/s) 

(4) 

Peak 

Velocity 

Shift 

(km/s) 

(5) 

Centroid 

Velocity 

Shift 

(km/s) 

(6) 

Velocity 

Dispersion 

(km/s) 

(V) 

Pearson 

Skewness 

Coeff. 

(8) 

Kurtosis 

Coeff. 

(9) 

Integr. 

Flux^ 

(10) 

EW 

(A) 

(11) 

001224 

2001/08/20 

3200 ± 300 

9300 ± 400 

-1880 ±90 

900 ± 200 

4200 ± 200 

0.23 ± 0.02 

2.86 ± 0.04 

3500 ± 100 

92 ±4 


2009/12/16 

2600 ± 300 

8700 ± 400 

-2080 ± 80 

6 ±200 

3900 ± 200 

0.26 ± 0.02 

3.13 ±0.04 

3200 ± 100 

104 ±4 


2010/10/12 

4400 ± 200 

8100 ± 300 

-2120 ± 60 

200 ± 200 

3900 ± 200 

0.19 ±0.01 

3.71 ±0.04 

9000 ± 200 

98 ±4 


2014/08/28 

4600 ± 300 

9000 ± 400 

-2030 ± 90 

400 ± 200 

4200 ± 200 

0.20 ± 0.02 

3.52 ± 0.04 

8300 ± 300 

134 ±6 

002444 

2000/12/22 

8490 ± 100 

13200 ± 200 

-170 ± 40 

100 ± 90 

4230 ± 70 

0.030 ± 0.008 

3.01 ± 0.03 

4010 ± 80 

93 ±3 


2009/12/18 

9400 ± 200 

14800 ± 400 

-410 ± 80 

-400 ± 200 

4800 ± 200 

0.01 ± 0.01 

3.14 ±0.04 

1840 ±50 

106 ±4 


2014/08/29 

10100 ± 200 

17400 ± 300 

-330 ± 60 

1900 ± 100 

6600 ± 100 

0.155 ±0.009 

3.64 ± 0.03 

5400 ± 100 

136 ±5 

015530 

2001/09/16 

7500 ± 200 

10500 ± 300 

2120 ± 70 

1100 ± 200 

3500 ± 200 

0.04 ± 0.01 

3.43 ± 0.04 

6300 ± 200 

108 ±4 


2009/12/17 

7200 ± 200 

10100 ± 300 

2200 ± 70 

1000 ± 200 

3400 ± 200 

0.07 ± 0.01 

3.32 ± 0.04 

2460 ± 70 

98 ±4 


2011/08/31 

7100 ±200 

10000 ± 300 

2320 ± 60 

600 ± 200 

3000 ± 200 

-0.01 ± 0.01 

2.86 ± 0.04 

4300 ± 100 

77 ± 3 


2014/08/28 

6400 ± 300 

11500 ±500 

40 ± 100 

800 ± 200 

4400 ± 300 

0.06 ± 0.05 

3.55 ± 0.05 

7100 ±300 

103 ±5 


^ Integrated fluxes are in units of 10 erg s ^ cm" 


Table 4 

Spectral measurements for [O III] A5007 [Abridged] 


Object 

Name 

SDSSJ 

(1) 

Obsei'vation 
Date (UT) 
(2) 

A(5100A)“ 

(3) 

FWHM*' 
(km s“') 

(4) 

Velocity 
Dispersion 
(km s“^) 

(5) 

Pearson 

Skewness 

Coefficient 

(6) 

Kurtosis 

Coefficient 

(V) 

Integrated 

Flux" 

(8) 

EW 

(A) 

(9) 

001224 

2001/08/20 

36.0 ± 0.2 

620 

310 ±20 

-0.20 ± 0.04 

3.17 ±0.05 

360 ± 20 

10.0 ±0.6 


2009/12/16 

29.0 ± 0.2 

750 

510 ±20 

-0.14 ±0.04 

3.44 ± 0.05 

460 ± 20 

16.0 ±0.8 


2010/10/12 

84.0 ± 0.2 

680 

570 ± 10 

-0.07 ± 0.03 

3.80 ± 0.04 

1590 ± 50 

18.0 ±0.7 


2014/08/28 

59.0 ± 0.2 

780 

780 ± 20 

0.08 ± 0.04 

3.83 ± 0.05 

1360 ± 70 

23 ± 1 

002444 

2000/12/22 

39.0 ± 0.2 

410 

298 ±8 

-0.06 ± 0.02 

4.12 ±0.02 

544 ±9 

13.0 ±0.3 


2009/12/18 

16.0 ± 0.2 

520 

390 ± 20 

-0.05 ± 0.03 

4.82 ± 0.04 

260 ± 10 

16.0 ±0.8 


2014/08/29 

37.0 ±0.2 

490 

340 ± 10 

0.18 ±0.03 

4.53 ± 0.03 

510 ± 10 

13.0 ±0.4 

015530 

2001/09/16 

54.0 ± 0.2 

280 

160 ± 10 

-0.06 ± 0.03 

3.69 ± 0.04 

1460 ± 60 

26 ± 1 


2009/12/17 

23.0 ±0.2 

470 

250 ± 10 

-0.28 ± 0.03 

3.75 ± 0.04 

870 ± 30 

36 ± 1 


2011/08/31 

52.0 ±0.2 

310 

230 ± 10 

-0.16 ±0.03 

4.40 ± 0.04 

1310 ±40 

25 ± 1 


2014/08/28 

65.0 ±0.2 

470 

200 ± 30 

0.08 ± 0.04 

3.19 ±0.06 

1800 ± 100 

28 ±2 


^ Flux densities are in units of 10 erg s ^ cm ^ A f 
^ Uncertainties on the [O III] FWHM are less than 70 km s“^ 
Integrated fluxes are in units of 10^^^ erg s^^ cm^^. 


For either definition of the structure function, the variabil¬ 
ity can only be characterized if it is detected above the un¬ 
certainties on the measurements. In practice, this has the fol¬ 
lowing consequences. First, on short time scales where the 
variability amplitude is small, we are not always able to mea¬ 
sure the structure function. Second, because the two defini¬ 
tions of the structure function are slightly different, the tem¬ 
poral bins where the structure function is undetermined can 
differ. This means that for some time intervals, one of the two 
definitions yields a meaningful measurement of the structure 
function while the other does not. 

For Gaussian distributions of Am, SF and SFjqr will be 
identical, but for distributions with more prominent wings 
(e.g., an exponential distribution), the value of SFiqr will be 
smaller than the value of SF corresponding to the same time 
interval. SF/qr is particularly robust against outliers that re¬ 
sult from bad data and is an especially sensitive measure of 
variability at short time scales where the amplitude of the 
variability is small (iMacLeod et alJl2012h . Because we are 
particularly interested in longer time scales with the SBHB 


candidates, it is not clear that the IQR structure function is 
preferable, so we consider both forms of the structure func¬ 
tion. Given the distribution of magnitude differences in Fig- 
ures|2]and[3 we expect to find SFigniAt) < SF(At). 

In calculating the structure function, it is common to in¬ 
clude a term to account for the photometric uncertainty. This 
can be done by calculating the structure function either of 
nearby non-varying stars, or [Olll], which is assumed to re¬ 
main constant. However, we are assuming here that the [O III] 
flux remains constant and we use it to re-normalize the spectra 
before calculating Am, which, in turn does not allow us to de¬ 
termine the photometric uncertainty in this fashion. Instead, 
we propagate the uncertainties in the integrated [Olll] flux 
resulting from the spectral decomposition into the error bars 
on the structure function. In some cases these uncertainties 
can be very large, so we calculate the structure function using 
only measurements with small uncertainties (as determined 
from three iterations of 3 (t clipping on the distribution of per¬ 
cent uncertainties in each At bin). In practice, this approach 
trims a small number of data points with uncertainties greater 
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Figure 10. A comparison of the continuum (black, filled circles) and H/3 
(blue, open circles) structure functions for the SBHB candidates. The two 
versions of the structure function presented in Equations m and are 
shown in the top and bottom panels, respectively. At the shortest time scale 
the value of the structure function here is driven primai'ily by observations 
of one object and may not accurately reflect the ensemble variability of the 
sample. We find that the shapes of the continuum and emission line structure 
functions are very similar, i.e. the majority of data points match within the 
uncertainties. 

than 50%, leaving the majority of points with uncertainties of 
less than 10%. 

The continuum and H/3 structure functions for the SBHB 
candidates are shown in Figure (TO] We have selected loga¬ 
rithmic time bins so as to get an equal number of pairs of ob¬ 
servations in each bin and adopt this binning for all the struc¬ 
ture functions of the SBHB candidates in this work. We find 
that the structure functions for the continuum flux density and 
the integrated H/3 line flux agree with each other within the 
uncertainties for the SBHB candidates. This is true for both 
definitions of the structure function that we have adopted here. 

For both definitions of the structure function, the variability 
in the shortest time scale bin (< 20 days) appears very large. 
While most observations in this bin show small variability, 
the value of the structure function here is driven primarily by 
many observations of one object which shows large variabil¬ 
ity. Thus, we do not regard this point as indicative of the 
ensemble variability of the SBHB candidates and exclude it 
from later figures. 

3.3. Comparison to the Sample of Typical SDSS Quasars 

To assess whether the behavior of the SBHB candidates dif¬ 
fers significantly from that of typical quasars, we compare 
their variability properties to those of the sample of typical 
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Figure 11. A comparison of the continuum (black, filled circles) and H/3 
(blue, open circles) stru cture functions for the comparison sample. The 
IVanden Berk et all 1200411 and IQR structure functions are shown on the top 
and bottom, respectively. The shape of the structure function for the contin¬ 
uum and emission line is similar, but not always consistent with each other 
within the uncertainties. 

quasars presented in Section 122] We measure the structure 
functions of this sample in the same fashion as the SBHB 
candidates, although we use a new set of temporal bins for 
this purpose, appropriate for the time sampling of the quasar 
fluxes. 

The magnitude differences of the continuum and the H/3 
are shown as a function of the rest-frame time interval in Fig¬ 
ures |5] and |6] respectively. The time intervals between ob¬ 
servations are clearly concentrated at values below 700 days, 
much shorter than those available for the SBHB candidates. 
This is an inevitable consequence of the fact that the spectra of 
the comparison sample were obtained as part of the SDSS I- 
III surveys. In contrast, the follow-up spectra of many of 
our SBHB candidates were obtained long after their origi¬ 
nal SDSS spectra, therefore the time interval between those 
observations are longer. A cursory inspection of these figures 
shows relatively few objects exhibiting relative flux variations 
greater than 50%. The distribution of continuum magnitude 
differences, shown in the right panel of Figure|5] has a median 
of 0.03 and standard deviation of 0.22. The distribution of H/3 
magnitude differences, shown in the right panel of Figure [b] 
has a median of 0.02 a standard deviation of 0.28. We char¬ 
acterize the ensemble variability properties of the compari¬ 
son sampl e by calculating the structure function, as discussed 
in Section 13.21 The structure functions for the continuum at 
5100 A and the H/3 emission line are presented in Figure [TT] 
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Figure 12. A comparison of the 5100 A continuum structure functions the 
SBHB candidates (black, filled circles) and the comparison sample (red, open 
circles). The dashed line shows the modified exponential parameterization of 
combined comparison sample and SBHB structure functions. See Section |3] 
for additional details. 

According to the structure functions, the fractional variability 
amplitude in this sample is generally on the order of 10-30%, 
although it increases with increasing time interval. 

In Figure [12] we compare the continuum structure function 
of the SBHB candidates to that of the comparison sample and 
in Figure [13] we show an analogous comparison for the H/? 
structure functions. We find that the SBHB candidates dis¬ 
play variability that is consistent with that observed for typi¬ 
cal quasars in a similar redshift and luminosity range, on the 
order of 10-30%. On the shortest time scales the variability 
may be even less, but that measurement is hindered by the 
magnitude of the uncertainties. 

To ascertain the robustness of the structure functions and 
verify the error bars, we simulate the effect on the structure 
function of perturbing the data within the measured uncertain¬ 
ties. The Am distribution is redrawn by randomly selecting 
flux ratios for each object in the sample from a normal dis¬ 
tribution with mean and standard deviation equal to the flux 
ratio measurement and uncertainty, respectively. The Am dis¬ 
tribution is then propagated into a structure function assum¬ 
ing the observed uncertainties. We repeat this process 10^ 
times and determine the median of the SF and SFiqr distribu¬ 
tions in each time bin and their associated standard deviations. 
We And that this method reproduces the observed flux ratio 
distribution within the uncertainties and preserves the mean. 
We conclude that our uncertainties are reasonably character¬ 
ized and that the match between the structure functions of the 



10 20 50 100 200 500 1000 2000 
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Rest-frame time interval [days] 

Figure 13. A comparison of the H/3 structure functions for the SBHB can¬ 
didates (black, filled circles) and the comparison sample (red, open circles). 
The dashed line shows the modified exponential parameterization of com¬ 
bined comparison sample and SBHB structure functions. See Section[3]for 
additional details. 

SBHB candidates and the comparison sample is robust. 

We also ask what the structure functions would look like 
with zero variability, given the observed uncertainties. This 
is simulated in an identical fashion, except the flux ratios are 
generated from a normal distribution around unity. We find 
that zero variability gives a flat structure function near but not 
equal to zero. This occurs because the structure function in¬ 
volves the root-mean-square of the absolute values of magni¬ 
tude differences, therefore, positive and negative excursions 
from the mean magnitude do not cancel out. This confirms 
that, within the observed errors, the observed structure func¬ 
tions are not consistent with zero variability. 

In view of the different distributions of time intervals for 
the SBHB candidate sample and the comparison sample, we 
can only make a direct comparison of the measured structure 
functions for time intervals between 300 and 900 days. To 
assess the similarities between the structure functions of the 
two samples further, we ask whether the models in the litera¬ 
ture that are used to describe the structure functions of typical 
quasars can also describe the structure functions of the SBHB 
sample. 

The structure functions of typical quasars are generally 
well parameterized by s imple power-law fun ctions (e.g., 
IVanden Berk et al] l2004t iMacLeod et al.l 12012h as well as 
modified exponential functions that are motivated by damped 
random walk models. Because we see evidence of a turn-over 
in the structure functions of the SBHB candidates on time 
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Table 5 

Structure Function Parameterization^ 


Type 

(1) 

SToo 

(mag) 

(2) 

T 

(days) 

(3) 

X^lv 

(4) 

Continuum 

0.324±0.026 

360±150 

2.8 

Continuum IQR 

0.289T0.017 

410±90 

5.9 

H/3 

0.335T0.021 

8±6’’ 

2.1 

H/3 IQR 

0.429T0.058 

1560T590 

8.1 


“ The modified exponential function is defined in 
Equation |4] All parameterizations listed here are 
for the combined binary candidate and SDSS com¬ 
parison samples. 

T his va lue is not constrained, as discussed in Sec¬ 
tion |33] 

scales around 1000 days, we adopt a modified exponential 
here. There is some variation in the exac t form that is adopted, 
for example IVanden Berk et all (120041) adopt a power-law 
model, but we choose the form used bv iMacLeod et al.l (12012h 
because it is able to reproduce the break that we see in our 
structure functions. The parameterization is 

1 /2 

5F = 5Foo(l-e-^'/") , (4) 

where SFoo and r are adjustable parameters describing the 
normalization and break time, respectively. The resulting fits 
are also shown in Figures [12] and [13] with the best-fit param¬ 
eters listed in Table [5] Note that we find r ^ 10 days for 
the H/3 structure function calculated from Equation]!] (shown 
in the top panel of Figure [T3ll. This value is not constrained 
because the variability on short time scales cannot be distin¬ 
guished above the noise, hence no turnover is observed in the 
structure function. 

We find that the parameterizations of the combined SBHB 
candidate and comparison structure functions here are gener¬ 
ally consistent with those in the literature. For structure func¬ 
tions of the continuum light between 4000 and 6000 A, sam- 
pled on rest-frame time scales of approximately 1-1000 days, 
livtacFeod et al.l (120121) find that their data are well fit by 
SFac = 0.222 ±0.13 mag and r = 638 ± 98 days. This break 
and normalization are comparable to the parameterizations for 
both continuum structure functions. 

In the range of time intervals from 300 to 900 days, the 
structure functions of the SBHB candidates and the compar¬ 
ison sample overlap and are in good agreement with each 
other. Furthermore, the structure functions of the SBHB can¬ 
didates and comparison sample can be well fit with a common 
model described by parameters comparable to those obtained 
for samples of normal quasars. This suggests that the SBHB 
candidates exhibit variability consistent with that expected for 
regular quasars on timescales > 900 days. Hence, we con¬ 
clude that the variability properties of the binary properties 
are similar to those of regular quasars. 

4. IMPFICATIONS OF OBSERVED VARIABIFITY 
PROPERTIES 

Our main, new observational result is that, on time scales 
longer than approximately one year the variability of both the 
continuum flux and the broad H/3 emission line flux in the 
SBHB candidates is consistent with the 10-30% observed in 
typical quasars of similar redshift and luminosity. Addition¬ 
ally, the continuum structure functions of the SBHB candi¬ 


dates and the quasars of the comparison sample can be de¬ 
scribed by a common model. The same is true of the H/3 struc¬ 
ture functions. Since we are only able to sample the stochas¬ 
tic variability of the SBHB candidates in a statistical sense 
and on relatively short time scales, we cannot carry out direct 
tests of specific models for the int eraction of a SBHB w ith a 
circumbinary disk . The models of iD’Orazio et al.l (120131) and 
iFarris et aH (120141) predict periodic variability on time scales 
ranging from about half to about 10 times the orbital period. 
The variability arises from accretion rate fluctuations on the 
orbital period which can also lead to flux variability because 
of beaming with half that period. On time scales longer than 
the orbital period, the variability is caused by tidal interactions 
of the binary with the circumbinary disk, leading to bright 
spots in the disk. These time scales would amount to several 
decades to many centuries for the SBHB candidates in our 
sample, in sharp contrast to the time scale we have probed, 
which are ^ 1 month to ^ 12 years. In order to carry out direct 
tests of the models, we would need to obtain light curves of 
individual objects spanning much longer temporal baselines. 

These results suggest three possible conclusions, (a) The 
binary candidates are, in fact, not SBHBs, and the mechanism 
responsible for the offset broad lines does not leave an imprint 
on the variability properties of these objects. In this case, the 
binary candidates still represent a very interesting population 
of quasars where some mechanism in the BLR has produced 
offset broad emission lines, as well as the correlation noted in 
Paper I and reproduced here where the lines are skewed in the 
direction of their shift. For example, such a mechanism may 
be a 1-arm spiral or a warp that leads to a high-contrast, non- 
ax isymmetric perturbation in a disk -like BLR (see examples 
in lLewis et al.lhoi^lWu et al.ll2008l respectively). Thus, the 
peak of the broad line is displaced according to the projected 
velocity of the brightest portions of the disk, (b) Some objects 
in the sample are SBHBs, but the contamination of our sam¬ 
ple by objects that are not is enough to prevent our detecting 
a difference, even though SBHBs have different variability 
properties, (c) SBHBs have the same variability properties 
as regular quasars. In order to say more about the first two 
points, we require the accumulation of additional lines of ev¬ 
idence, particularly the radial velocity curves that will be the 
focus of an upcoming paper. 

In the context of the third possibility, our main result sug¬ 
gests that, if the observed H/3 flux variability is the result of 
reverberation of the continuum, then the extent of the BLR of 
the SBHB candidates is comparable to that of typical quasars. 
Thus, in the remainder of this section, we consider the im¬ 
plications of these results for the structure of the BLR in the 
SBHB candidates. We restrict our discussion to the specific 
scenario in which the offset broad emission lines arise in a 
BLR that is asso ciate d with the sec ondary BH in a SBHB. 
IShen et aTl (1201 3h and lTu et^ (1201 3h have already attempted 
to constrain the orbital properties of the SBHB population 
based on the observed properties of the candidates. Here, we 
are interested in examining a more specific question in the 
context of our variability data: if the BLRs of SBHB candi¬ 
dates have a similar extent to those of typical quasars, what 
are the ensuing constraints on the orbital parameters of the 
hypothesized SBHBs? To this end, we compare the extent of 
a typical BLR to characteristic length scales in the binary, the 
size of the Roche lobe of the accreting black hole and the ra¬ 
dius at the which the accretion disk of the accreting black hole 
is truncated because of the tidal effect of its companion. 

We begin by collecting here a number of constraints derived 
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by combining observations and basic orbital mechanics. As 
in Paper I, we assume that only the secondary (less massive) 
BH is active, in other words, all observed properties of the 
AGN, such as the bolometric luminosity and the Eddington 
ratio, refer to the secondary. Adapting equations (1) and (2) 
of Paper I, we write the orbital period and separation of the 
binary as 

. 2625 Mg 

p = - 5yr (5) 


and 


0.432 Mg 


pc, 


( 6 ) 


where q = M 2 /Mi < 1 is the mass ratio. Mg is the total mass 
of the binary in units of 10* M©, and 14,3 is the true orbital 
velocity of the secondary in units of 10* km s“', assuming a 
circular orbit. In Paper 1 we presented the observed distribu¬ 
tion of projected orbital velocities, M 2 = V 2 sin/sin^ (where 
/ and (j) are the inclination and phase angles of the binary), 
spanning a range from one to a few thousand km s“'. There¬ 
fore, as long as / and do not have extreme values, we ejmect 
14,3 to be of order a few. Combining equations (|5]l and (|6]) we 
can cast Kepler’s third law for such a binary as 



/ p \2/3 

a = 0.228Mf pc. (7) 

In Figure [14] we show graphically the physical parameter 
space described by the binary separation, a (in pc), and the 
total mass of the two BHs, Mg(= M,o,/10* Mq). The dot¬ 
ted lines in Figure [T4| show the relation between a and Mg 
(a cx Mg; see eqn. |6]l, for a fixed value of the orbital speed 
of the secondary BH (V 2 = 4000 km s“*), and three different 
values of the mass ratio {q = 0.01, 0.1, and 1). These lines 
effectively delineate the entire range of combinations of a and 
Mg since the q = 0.01 line is nearly identical to the q = Q line 
and since q > \ contradicts our basic assumption that only 
the secondary BH is active. The solid lines in Figure [T4l are 
lines of constant orbital period, following Kepler’s third law 
{a cx Mg' ; see eqn|7]l. If the orbital speed of the secondary BH 
increases, then the lines of constant q move towards the lower 
right; their horizontal shift corresponds to a change in Mg by 
a factor (V 2 / 4 OOO km s“*)* as prescribed by equation (|6]l. 

In order for gas in the vicinity of the secondary black hole to 
be gravitationally bound to it, the gas must be within its Roche 
lobe. The following parametric expression for the effective 
radius of the Roche lobe of the secondary, R 12 (the radius of 
a sphere with the same volume as t he Roche lobe), r elative to 
the binary separation, a, is given bv lEgglet^ (119831) 


0.49 g*/* 

a Q.6q^P+ \n{\+q^P) 


( 8 ) 


At some critical distance from the secondary, the accreting 
object stops dominating the potential, which departs signif¬ 
icantly from spherical symmetry with the result that orbits 
(coplanar with the binary orbit) at that radius cannot be cir¬ 
cular. Therefore an accretion disk around the secondary is 
truncated at this critical distance because particle orbits start 
to cross each other. Detaile d calculations of the disk trunca¬ 
tion radius are presented by iPaczvnsldl (119771 and references 
therein). For our purposes we adopt the following approxi¬ 
mate expression for the truncation radius relative to the Roche 


Figure 14. Constraints on the binary parameters of the SMBH candidates 
derived by considering the size of the BLR in compaiison to other character¬ 
istic lengths in the binary. The dotted lines are lines of constant mass ratio for 
a fixed orbital velocity of the secondary BH. The thin solid lines are lines of 
constant orbital period. The conditions that the BLR is not tidally truncated 
and is smaller than the Roche lobe are depicted by grey bands. More details 
can be found in Section|4|of the text. 


lobe radius bv lEgglet^ (1201 IL section 6.1, eqn 6.8) 

^ ^ (l.9+—')4(^), (9) 

a \ q ) 

where xijq) is given by equation (|8]l. In Figure [15] we illus¬ 
trate how the secondary Roche lobe radius and the truncation 
radius vary with the mass ratio for a binary with a = 0.3 pc 
and Mg = 1 (hence P sa 1500 yr). For a larger value of a, R 12 
and R,r increase proportionally. However, R,r is always small 
compared to R 12 and a. 

To estimate the extent of the BFR we start wit h the empir¬ 
ical ra dius-luminosity (Rblr-L) relationship of iBentz et al.l 
(1201 3h . connecting the reverberation lag of the broad H/3 flux, 
th/ 3 , to the monochromatic luminosity of the AGN at 5100 A. 
This relationship is calibrated over four orders of magnitude 
in monochromatic luminosity and has some scatter, whereas 
our sample spans a relatively narrow luminosity range (about 
an order of magnitude). To account for the scatter, we do 
not apply the Rblr-L directl y but instead we c onsider the 
measured lags of the AGNs in iBentz et al.l (1201 3h with lumi¬ 
nosities in the same range as the AGNs of our sample. For 
this narrow range of luminosities there is no discernible cor¬ 
relation therefore we simply adopt the median lag to obtain 
an estimate of Rblr = cth/? = O.On/j^ QQg pc, corresponding to 
a monochromatic luminosity of log[ALA(5100 A)/ergs“'] = 
43.62. The upper and lower bounds on Rbir define the in- 
tervals within which we find 1 /3 of the iBentz et aH (120131) 
AGNs above and below the best value of Rblr, respectively 
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(analogous to “ilcr” intervals). The extent of the BLR (its 
effective outer radius) is likely a few times larger than Rblr, 
as evidenced by the following observational results: (a) The 
Ha line responds to continuum variations with a lag that is 
approximately twice as long as that of Wf3 (e.g.. iGrier et al.l 
120131 iBentz et al.ll2010h . Similarly, the lines of the optical 
Fe II complex respond to continuum variations wi th a lag that 
is app roximately 2-3 times as long as that of H/? (iBarth et al.l 
120131) . (b) The dust sublimation radius, as inferred from inter¬ 
ferometric imaging and from the reverberation lag of the near- 
IR continuum relative to the optical continuum is 4-5 times 


larger than Rbtr (e.g 


limoto et alJl201 ItlSuffanuma et al.l 

lIoMiKoshida et al .1120 14 

). In view of the above results, we 


take Rdust = ^Rblr to be a plausible estimate of the maximum 
extent of the line-emitting gas in the BLR of a typical quasar. 
Because of the scatter in the Rblr-L relation we illustrate the 
outer extent of the BLR as a broad band rather than a sharp 
line in Figure [14] In the illustration of Figure [15] we show 
the BLR and dust sublimation radii as dashed lines. For the 
combination of binary parameters adopted in that illustration, 
Rblr can be either smaller or larger than R,r depending on the 
value of q, but it is always smaller than Ri 2 - In contrast, Rdust 
is always larger than Rtr, but it can be either smaller or larger 
than Ri 2 depending on the value of q. 

If we suppose that the extent of the BLR of SBHB candi¬ 
dates is the same as that of typical quasars, then we obtain 
lower limits on the orbital separation and total mass in the bi¬ 
nary. We consider two scenarios to place constraints on the 
binary properties, amounting to two ways of placing the trun¬ 
cation radius relative to the dust sublimation and BLR radii 
(see example in Figure [TsTl. 

In the first, more extreme scenario, we require that 
Rdust ^ Rtr, which is equivalent to assuming the BLR has the 
same extent as in a typical quasar of a given luminosity. This 
condition can also be viewed as a requirement that the po¬ 
tential in the outskirts of the BLR is still dominated by the 
secondary BH. If we adopt this scenario, we effectively as¬ 
sume that the binary separation is large enough that the BLR 
must not be tidally truncated by the primary black hole. Thus 
we obtain the dashed black line at the upper right corner of 
Figure [14] which is based on the best estimate of Rdust- The 
grey band represents the uncertainty in Rdust about the best 
estimate, arising from the scatter in the correlation between 
Rdust and quasar luminosity. The allowed values of a and Mg 
are in the strip between the extreme values of q and above the 
Rtr = Rdust band. This leads to larger orbital separations and 
we infer large total masses, 10'° Mq, and very long orbital 
periods, ^10^ yr. 

In the second scenario we require that Rdust ^Rl 2 , which is 
equivalent to assuming only that the BLR is bound to the sec¬ 
ondary BH. This case is more conservative in the sense that 
it results in less stringent constrains on the properties of the 
binary given that the extent of the BLR is set by the observed 
quasar luminosities and R12 > Rtr- This more generous as¬ 
sumption leads to masses .^10® Mq and periods .^10^ yr, as 
indicated by the dashed black line and grey band at the lower 
left corner of Figure [14] We consider this possibility more 
likely because the emission line profiles of our SBHB candi¬ 
dates are generally asymmetric, unlike the typical quasar pop¬ 
ulation (see, for example. Figure 5 of Paper I). Their asym¬ 
metries suggest that the structure of the BLR in SBHB can¬ 
didates is not typical; this may result from the tidal effect of 
the primary BH onto gas that is bound to the secondary BH. 
The above constraints are modified if we adopt a higher value 



Figure 15. Comparison of the secondary Roche lobe radius (R/, 2 . solid black 
line) and truncation radius {Rtr, solid red line) to the empirically calibrated 
BLR radius (Rblr, lower dashed line) and dust sublimation radius (Rjust, up¬ 
per dashed line). In this particular realization we have assumed a binary with 
a total mass of 10* M 0 (Mg = 1) and an orbital separation of a = 0.3 pc and we 
plot Rl 2 and Rtr as a function of the mass ratio, q. We take Rblr = 0.017 pc 
and Rdust = 0.068 pc; these are fixed since they depend on the luminosity of 
the active BH. More details can be found in Sectionjdjof the text. 


of the orbital speed of the accreting BH. Specifically, we can 
set V 2 = 8000 km s“* without violating the observational con¬ 
straints, which would lead to an allowed mass range that is 
four times higher than that depicted in Figure [14] for a given 
orbital separation, and orbital periods that are approximately 
a factor of two shorter. 

The question of the allowed orbital periods is important be¬ 
cause it has a bearing on the detectability of radial velocity 
variations. If the expected orbital periods are of order 10^ yr, 
as in the first scenario described above, then we would not 
be able to detect significant radial velocity variations in our 
ongoing spectroscopic monitoring program with observations 
spann ing a decade (see illustration in Fig. 5 of iDecarli et al.l 
l2013h . Thus, if we do detect significant velocity variations, 
we may conclude that, in the context of the SBHB scenario 
one or more of the following may be true: (a) the mass of 
SBHB is ^several x 10^ Mq, as long as the secondary BH 
is the one that is active, (b) the BLR is not confined within 
the truncation radius defined by equation (]9]), (c) the BLR is 
associated with the primary BH. Possibility (c) implies that 
^ > 1, effectively, which corresponds to the lower right cor¬ 
ner of Figure [T^ where the equations above no longer apply 
and must be re-cast accordingly. 

5. SUMMARY 

We have presented new spectroscopic observations of a 
sample of 88 quasars that are candidate SBHBs. These were 
selected in Paper I based on the fact that their broad H/3 emis¬ 
sion lines are displaced from the frame of their host galaxies 
by rsu 1000-5000 km s“\ which was taken as the signature 
of motion of the accreting BH around an unseen compan- 
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ion. Moreover, in the first round of followup observations, 
presented in Paper I, ~ 20% of the objects observed showed 
changes in the velocity offsets of their broad H/3 lines that 
were consistent with binary motion. Our primary goal in this 
paper was to present the new observations and improvements 
in our analysis methodology, including a more extensive anal¬ 
ysis of the uncertainties. The new observations extend the 
temporal baseline of our ongoing monitoring campaign and 
form the basis for a study of radial velocity variations, which 
will be the subject of a forthcoming paper. 

We have also taken this opportunity to study the variability 
of the optical continuum and the integrated flux of the broad 
H/3 emission line in the sample of SBHB candidates. We 
compared the variability properties of this sample with those 
of a comparison sample of typical quasars of similar redshift 
and luminosity with multiple epochs of spectroscopy avail¬ 
able from the SDSS. We decomposed the optical spectra to 
separate the quasar continuum, broad H/3, and narrow [O III] 
emission and calculated spectral properties. These measure¬ 
ments, which are of general interest for the sample but are not 
all used in this work, are presented in an appendix along with 
a careful analysis of the associated uncertainties. 

For each sample, we calculated the structure function for 
the optical continuum and broad H/3 emission line. We found 
that the ensemble variability of the SBHB candidates is con¬ 
sistent with that of the comparison sample of typical quasars 
on time scales of approximately one year or longer. Moreover, 
the continuum and H/3 structure functions of the two samples 
are described by a common model, with parameters applica¬ 
ble to typical quasars. This result is robust in the sense that it 
cannot be a consequence of measurement uncertainties and is 
inconsistent with no variability. Furthermore, the amplitude 
of variability of the broad H/3 emission line on a given time 
scale is comparable to that of the optical continuum measured 
at 5100 A. One possible conclusion is that, taken together and 
assuming that the H/3 flux variability is driven by the con¬ 
tinuum variability, these results suggest that the extent of the 
BLR is SBHB candidates is similar to that in typical quasars. 
But other interpretations are possible, as we note at the begin¬ 
ning of section |4] 

We have considered the implications of these results in the 
context of the binary hypothesis. If we suppose that the BLR 
of SBHB candidates is similar to that of typical quasars and 
is not truncated by the binary companion (our first scenario), 
we infer long orbital periods (.J/, 10^ yr) and large total masses 
10'*’ Mq). If, however, we only require that the BLR gas 
is contained within the Roche lobe of the accreting BH and 
its outer parts can feel the tidal influence of the companion 
(our second scenario, a more likely possibility in view of the 
skewed and asymmetric H/3 profiles), we infer orbital peri¬ 
ods and masses that are approximately an order of magnitude 
lower. These inferences can be tested by measurements or 
limits on radial velocity variations, which will be the focus of 
a forthcoming paper. 
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APPENDIX 

EVALUATION OE UNCERTAINTIES IN SPECTRAL MEASUREMENTS 

In this appendix we carry out a variety of tests to quantify the uncertainties in the quantities reported in this paper, including the 
moments of the line profiles and other properties of the broad emission lines. We pay particular attention to uncertainties arising 
from the decom position of the observed spectra since the method one uses for the decomposition may lead to systematic errors in 
the results (e.g.. lDennev et allUOllURunnoe et al.ll2013l) . We begin by comparing the results from two different decomposition 
methods, the one adopted in this paper and the one adopted in Paper I. We then consider how the measurement of moments of 
line profiles is influenced by noise in the wings of the profiles, which leads us to define an optimal spectral window around the 
peak of the line profile for making such measurements. Einally we examine the differences resulting from making measurements 
on the observed spectra and on the parametric models that describe these spectra and close with a comparison of the profiles 
moments measured in Paper I and those measured here. 

The Effect of the Spectral Decomposition 

iDennev et al.l (l2009h considered this issue for two objects observed multiple times during reverberation mapping campaigns 
and concluded in those cases that the measurements yield similar results. We are interested in this question generally, but also 
more specifically for objects with shifted broad emission lines. 

We use the SDSS spectra from our sample to address this question because we have d ecom posed them twice, first following the 
method described in Paper I and again following the new method describe d in Section lZ3] There are two significant differences 
between these methods. Eirstly, Paper I uses the lBoroson & GreenI (Il992h Ee II te mplate whereas here we ad opt the template of 
IVeron-Cettv et al.l ([2003). The primary difference between these templates is that IVeron-Cettv et al.l (120041) identified a region 
of narrow lines and excluded these from their Ee ll template. In practice, the two templates are similar with a small difference 
around 5000 A that can be important, particularly when the Ee ll emission lines are strong and/or narrow. Secondly, we adopt a 
different procedure for subtracting the narrow lines than in Paper I. In both cases, an [O III] line is used to constrain the shape of 
the narrow In Paper I, a narrow line template was constructed by fitting the [O Hi] profiles with a cubic spline after fitting the 
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Figure 16. A compari son o f two spectral properties measured in identical fashion on broad H/3 spectra that were decomposed following the two different methods 
described in Appendix lA. II Top: A direct comparison of the measured properties. The dashed line shows where the measurements are equal and red open p oints 
are rejected from statistics in Table [^by up to 5 iterations of 3a clipping. Middle: distribution of the measurement discrepancy described in Appendix EU] 
Bottom: The discrepancy as a function of S/N in broad H/3. The dashed line shows the median value for the discrepancy after clipping. 


underlying wing of the broad H/3 line with a low-order polynomial. Here we ht the [O III] prohles with combination of Gaussians 
instead. In practice, the method adopted in Paper I tends to subtract more of the underlying flux in the [O III] region than the one 
adopted here because no assumption is made about the shape the [O III] wings. 

In order to make measurements of the spectral quantities of interest (line widths, shifts, fluxes, and profile models; see Table|6l), 
we isolate the broad H/? profile by subtracting the continuum, and the Fe II, [O III] and narrow H/3 lines. The measurements are 
made directly on the data following the procedure adopted in Paper I and also on the best-fitting model for the broad H/3 profile 
as in the procedure adopted in this paper. The resulting measurements for two representative quantities are compared Figure [Tbl 
In this figure, the top row shows the measurements from each decomposition method plotted against each other with a unit slope 
line superposed for reference. The middle row shows the distribution of discrepancies between the measurements. We adopt a 
convention where the discrepancy in a quantity X is always taken to be AX = Vp^per i - Vpaper ii and the fractional discrepancy is 
Ax/{X). The bottom row shows the relation between the discrepancy between measurements and the S/N at the peak of the 
broad H/3 line. Statistical measures describing these distributions are summarized in Table|6] 
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Table 6 

Measurement discrepancies between spectral decomposition methods “ 


Property 

(1) 

Mean 

(2) 

Median 

(3) 

Standard 

Deviation 

(4) 

Fractional 
Difference ^ 

(5) 

Outlier 
Fraction ^ 
(6) 

Centroid Velocity Shift (km s“') 

-300 

-200 

300 


0.02 

Velocity Dispersion (km s“^) 

5 

50 

200 

0.08 

0.05 

Skewness Coefficient 

-0.08 

-0.04 

0.18 


0.10 

Pearson Skewness Coefficient 

-0.02 

-0.01 

0.03 


0.10 

Kurtosis Coefficient 

0.10 

0.06 

0.20 

0.07 

0.11 

FWHM (km s^‘) 

-100 

-70 

500 

0.09 

0.09 

FWQM (km s-‘) 

-600 

-500 

1000 

0.13 

0.02 

Peak Velocity Shift (km s“*) 

-30 

-20 

300 


0.06 

Int. flux (erg s“^ cm“^) 

-60 

-40 

100 

0.11 

0.07 

EW(A) 

9 

8 

10 


0.03 

Fe II int. flux (erg s“^ cm“") 

70 

70 

200 


0.05 

Fe II EW (A) 

7 

6 

11 


0.02 

[0 III] int. flux (erg cm“^) 

-4 

-1 

40 


0.06 

[0 III] EW (A) 

5 

3 

7 


0.07 

/x(5100 A) (erg cm“^ A“’) 

-2.9 

-2.3 

2.7 

0.20 

0.06 


“ Statistics after outlier rejection for the distributions of measurement discrepancies 
shown in Figurefl^ The discrepancy in a quantity X is taken to be AX = Xpap i -Xp-^p n. 

The listed value is the standard deviation of the fractional difference distribution, 
where the fractional difference is AX/(X). 

The outlier fraction is the fraction of objects rejected after up to 5 iterations of 3 ct 
clipping. 

We find that the distributions of discrepancies between measurements often have a small number of significant outliers, which 
we reject before calculating statistics for the distributions. We iteratively reject > Str outliers relative to the median until one of 
the following three conditions is met: a) there are no > 3o- outliers left b) fewer than 1% of the total points have been rejected in 
the last iteration, or c) we have repeated the process hve times. The rejected points are shown as red open circles in the relevant 
scatter plots and are included in the relevant histograms. These points are excluded when calculating the statistics reported in 
Table|6]and their fractional number is included in the last column of the table. 

To investigate whether the magnitude of the discrepancies depends on the S/N the peak of the broad H/3 line, we evaluate the 
root-mean-square dispersion about the best-fitting parametric model near the peak of the line and plot it against the discrepancy 
in the bottom row of Figure [161 In most cases displayed in these figures, the discrepancies between methods are smallest at the 
highest S/N and become larger at lower S/N. 

The measurements based on different spectral decompositions agree reasonably well. Trends in the measurement discrepancies 
(i.e., asymmetries in the histograms of Figure [Thll can be attributed to the difference in the [Olll] subtraction. The method of 
Paper I depresses the wing of the broad H/3 line more than the method of Paper II and makes the profiles cuspier. Thus the 
measurements from Paper I give bluer centroid wavelengths, somewhat larger skewness coefficients, and smaller FWQM for the 
broad H/3 lines. Moreover, the [Olll] lines in Paper I have somewhat larger EWs. The relative velocities measured by cross- 
correlating broad \{(3 spectra are not affected, although the absolute velocity of the broad H/3 lines measured from the original 
SDSS spectra incurs an uncertainty of 300 km s“* as a result of the different methods of fitting the prohles around the peak in 
Papers I and H. 


Practices for Measuring Spectral Properties 

The practices adopted for measuring spectral properties can vary between authors. Most notably, some authors make the 
measurements directly from the observed spectra wh ile oth ers ht parametric m odels to the spectra and make the measurements 
from the models (e.g.. lVestergaard & Petersonll2006l versus IShang et alJ|2007h . Another variation in practices is the wavelength 
around a broad emission line prohle used to make the measurements. We determine the best wavelength window for making 
spectral measurements by systematically evaluating measured spectral properties for two windowing methods: a fixed window, 
where a hard wavelength limit is adopted for all objects, and a flexible window that extends from peak wavelength until the flux 
has dropped to some percentage of the line peak. For the model spectra, we And that the best window is, 4600-5100 A for a 
fixed-width scheme and down to 1% of the peak flux for a flexible scheme. We prefer the window based on a percentage of the 
peak flux because of its flexibility; the window can expand to accommodate very broad lines, asymmetric lines, or, in the case 
of this sample, very shifted lines. In the case of measurements made on the data, a more “conservative” range is required for 
the flexible window so that the measurement is not influenced by noise. When making measurements on the observed data, it 
is therefore best to use a window going out from the line center until the spectrum reaches 7% of the peak flux. Measurements 
involving the height the line peak (e.g., FWHM) are not particularly sensitive to the window, whereas the moments of the line 
profile and integrated flux do depend on the window; the consequences of using too wide a window are much more significant 
in measurements made on the data. Following the methodology adopted for comparing measurements from different spectral 
decompositions, we investigate the differences between making measurements on the data versus the best-fitting model. In 
general, though it appears that measurements made on the data are more sensitive to the characteristics of and spikes in the 
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Figure 17. A compaiison of the velocity shift of the peak and Peai'son ske\vness coefficient measured for the SDSS spectra in this work and in Paper I. 

noise, we find that the overall effect is negligible (e.g., the fractional difference in the velocity dispersion is 2%) compared to the 
differences that can result from different spectral decomposition methods. 

A Comparison of Measurements in This Sample 

In this work we measure spectral properties on the best-fitting models from the decomposition described in Section 12.31 in 
the wavelength window where the model is greater than 1% of the peak model flux. This differs from Paper I in the spectral 
decomposition, the wavelength window of interest, and the spectrum on which the measurements are made. As a result, there is 
scatter between the properties originally measured for the SDSS spectra and those we adopt here. Measured spectral properties 
are listed here in Tables [3] and |4] and can be compared to the original values in tables 3 and 4, respectively, in Paper I. For the 
velocity shift of the peak and the Pearson skewness coefficient, we compare the Paper I measurements with our adopted values 
in Figure [TTl These properties are of particular interest because the velocity shift of the peak anchors the radial velocity curves 
and because of the correlation between the Pearson skewness coefficient and velocity shift of the peak identified in Paper I. The 
correlation shown in Figure 5 of Paper I illustrates the tendency for profiles to be asymmetric in the direction of their shift. We 
find that this correlation is robust against the details of the spectral decomposition and measurement and show it here in FigurefTSl 
For the measurements in this work, the Spearman Rank correlation coefficient is p = 0.419, with a probability P = 1.1 x 10“^ that 
this distribution of points would be found by chance. Thus, while the details of the correlation do depend on the measurements, 
the general trend is robust and statistically significant. 



Figure 18. The correlation between Pearson skewness coefficient and the velocity shift of the peak found in Paper L Open red circles show the measurements 
from Paper I and the solid black circles show the re-measured values from this work. 







